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Abstract 

In the presence of a low scale seesaw of type I + III, flavor violating effects 
in the leptonic sector are expected. Their presence in the charged sector 
is due to the mixing of the fermionic vector-like weak triplets with the 
chiral doublets, which cause non-universality of the tree-level Z coupling. 
We investigate the bounds on the Yukawa couplings which are responsible 
for the mixing and present the results for two minimal cases, a fermionic 
triplet with a singlet or two fermionic triplets. Different channels for these 
processes are considered and their current and future potential to probe 
these couplings is discussed. 



I Introduction 

Available experimental data on neutrino oscillations indicates a small mass 
of left-handed neutrinos. This is in contrast with the Standard Model (SM) 
where the left-handed neutrinos are massless. Also, the nature of neutrinos, 
whether they are Dirac or Majorana particles, is not known. The latter 
possibility is theoretically most compelling, since it introduces new physics 
at the scale A, where the neutrino mass operator 

or 5 = L -^- (i) 

is formed. There are only three different ways to realize this operator at the 
tree level when a single representation is added pQ. 

Adding a right-handed neutrino is referred to as the type I seesaw [21 
IH El [6], while an extra bosonic triplet with hypercharge 1 results in type 

II seesaw [3 El |9]. The third option is to couple the leptonic and Higgs 
doublets to a fermionic weak triplet with zero hypercharge and this is the 
type III seesaw [TUj . 
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Unfortunately, the scale A is not known since it depends on the size of 
the Yukawa couplings. If they are of order one, as in certain Grand Unified 
Theories (GUTs), Eq.Q predicts A around 10 13 GeV. Such a high scale 
would make it very hard to probe the origin of the mass operator directly at 
a collider. On the other hand, when Yukawa couplings in Eq.Q are small, 
the seesaw scale may lie anywhere below 10 13 GeV. Notice also, that the 
small Yukawa couplings are technically natural due to a protective chiral 
symmetry. 

Recently, a grand unified model has been proposed [TTJ which is an exten- 
sion of the minimal Georgi-Glashow SU(5) model with a fermionic adjoint 
representation, that predicts a low mass for a fermionic triplet from unifica- 
tion requirements [12] . Neutrino masses are realized with a combination of 
type I and III seesaws together with an upper bound on the fermionic triplet 
around TeV. Due to the fact that they are coupled to electroweak gauge 
bosons, one can produce the fermionic triplets at the LHC and measure the 
origin of neutrino masses by studying their decays [121 US EE]- In principle, 
one can even distinguish various seesaw types at a collider by studying events 
by their charged lepton multiplicity [15] . 

Besides the neutrino mass operator, also higher dimensional operators 
are produced below the seesaw scale, the size of which can be constrained 
by flavor changing processes [16] . While the leptonic mixing matrix becomes 
non-unitary in both type I and III cases [17] , the unique feature of the type III 
is the presence of charged lepton flavor changing neutral currents (FCNCs) 
at the tree level. Our aim is to establish, whether testing such processes may 
shed some light on the origin of neutrino mass in a minimal model with a 
predicted light triplet below TeV. 

In contrast to previous studies [18] [19], we use the existing data from 
neutrino oscillation experiments to express Yukawa couplings which enter 
the expressions for flavor violating processes. In other words, we relate the 
d = 6 operators with d = 5 by using a convenient parametrization [2TJ . 
The natural values of the Yukawa couplings for a low seesaw scale around 
the electroweak scale ~ 100 GeV are of the order of 10~ 6 . However, it 
turns out that there exists a portion of parameter space where the effect 
in these processes is observable, while neutrino masses remain small due to 
cancellations. 

We investigate current bounds on these couplings from various processes 
in two minimal cases with two heavy neutrinos, a singlet and a triplet, and 
two triplets. It turns out that at least in these minimal cases, the number 
of parameters which specify the rate is reduced to a single real parameter, 
which is most constrained by the /x — e conversion experiments. We use this 
constraint to asses other possible channels and also comment on non-minimal 
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Parameter Best fit 3a 

Am^[10- 5 eV 2 ] 
lAm^lflO^eV 2 ] 
sin 2 9 12 
sin 2 9 2 3 
sin 2 6»i3 



7.65 


7.05- 


8.34 


2.40 


2.07- 


2.75 


0.304 


0.25- 


0.37 


0.50 


0.36- 


0.67 


0.01 


< 0.056 



Table 1: Parameter fits from oscillation experiments taken from [22J. 

models in the end. 

We start with a discussion of the type I + III seesaw model of neutrino 
masses in section [2j where we focus on the two minimal cases and discuss 
the Casas-Ibarra-Ross parametrization. In section |3j we constrain the free 
parameters of the model using a bound from fi — e conversion searches in nu- 
clei. Next, we consider a comprehensive list of other constraints in section [4], 
we comment on non-minimal models in section [5] and present our conclusions 
and an outlook on future experiments in section [6] The appendices contain 
a derivation of the couplings of light and heavy leptons in models with arbi- 
trary number of additional fermionic singlets and triplets and a calculation 
of f2 — > /17 process in such models. 

2 Minimal type I and III models 

Neutrino oscillation experiments can be explained by non-zero masses of 
the light neutrinos. The best fit of the mass squared difference and the 
mixing angles are given in Table [T] and constitute evidence for a nonzero 
mass. While the neutrino masses are bounded from above by beta decay 
searches and cosmology, the overall scale of the neutrino mass has not been 
established, therefore the lightest neutrino may still be massless. We will 
consider two minimal models which accommodate the oscillation data, one 
with a singlet and a triplet (motivated by a GUT) and the other with two 
triplets. The reason for this choice is minimality and predictivity. In both 
cases, the lightest neutrino is massless and there is only one Majorana phase 
which cannot be rotated away. 

The Lagrangian for a model with a fermionic singlet and a zero hyper- 
charge triplet can be written in the following way using the two component 
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Weyl spinors 

Ct = iL^D^U + ilfa^D^ + iT^D„T a + itfa^S 

- {y'/H !.,(', - y l s H T iT 2 US - y l T H T ir 2 T a T a Li^ + h.c. (2) 

- 1/2 (m T T a T a + m s SS) + h.c. 

where stands for the appropriate covariant derivative. After spontaneous 
symmetry breaking, we obtain the well-known seesaw formula for the light 
neutrino masses 

(mr = -t(yM + yJd). (3) 

K J 2 I m T m s I K J 

When only type III is considered, the second term is replaced by the Yukawa 
couplings and the Majorana mass of the second triplet. 

The same Yukawa couplings responsible for the d = 5 operator, also 
contribute to the d = 6 operators. For example, the presence of ys alters the 
couplings of the W to the neutrino, which means that Upmns is no longer 
unitary, while the yx mixes the charged leptons and therefore also affects 
the universality of Z boson couplings. The Feynman rules for the fermion 
couplings in presence of a singlet and a triplet are presented in appendix [A} 

In order to use the information from the oscillation experiments on the 
neutrino mass to reduce the number of parameters, we employ a useful 
parametrization [21] . For the two minimal cases, the neutrino masses are 
fixed, because the lightest neutrino is massless, therefore Am 2 determines 
the mass of the heaviest two. This parametrization specifies all the Yukawa 
couplings in terms of measurable neutrino quantities and a single complex 
parameter z for the case of normal (NH) 

y l T = —i\f2mr/v {l}i2\lm v 2 cos 2 + Ui^m^ sin z^j , (4) 
y l s = —i\l2msjv {^—U-a^/m^ sin z + Ui^yjm^ cosz^ , (5) 

and inverted (IH) hierarchy 

yfp = —i\/2m T /v (Un \fu£{ cos z + L^a/^I s i n > (6) 
y l s = —i\ / 2ms/v (^—Ui2\/m'{ sin z + t/j2A/ m 2 cosz ^j > (7) 
where U is the unitary PMNS matrix defined by the standard parametriza- 



4 



tion and (J) is the additional Majorana phase 

(C12C13 S12C13 s 13 e~ tl 

-S12C23 - Ci 2 S 2 3Sl 3 e i<5 C12C23 - Si 2 S 2 3Sl3e t5 S23C13 
■S12S23 — Ci3C23Si3e 4 ' 5 — C12S23 ~~ Sl2C23 s 13 eWS C23C13 

xdiag(l,e^,l). 

The size of the Yukawa couplings is determined by the complex z parameter 
and it increases exponentially with lm(z). In this case, the effects on the d = 
6 operators responsible for lepton flavor violating effects become visible, while 
at the same time neutrino masses remain small due to an exact cancellation. 
The higher the seesaw scale, the more severe fine-tuning is needed in order 
to produce a visible effect because the d = 6 operators scale as A" 2 while 
neutrino masses go as A -1 . 

The new particles may also be light, for example as predicted in [TT]. If 
this is so, we have the possibility to produce them at a high energy collider 
and measure the Yukawa couplings by decay rates and branching ratios [32] • 
In this paper we instead investigate various lepton flavor violating (LFV) 
processes and determine the values of lm(z) which are needed in order to 
observe them. 

When Yukawa couplings are large, e lra ^ factorizes, multiplies all the 
Yukawa couplings in Eqs.([6])-([7]) and there is no dependency on the real 
component of z. Therefore, we can state all the limits on d = 6 operators 
at a reference triplet mass with a single parameter, the imaginary part 
of z, which governs the overall size of d — 6. The ratios between the vari- 
ous channels are not affected by Im(z) and depend solely on neutrino mass 
parameters, together with the Majorana phase. This means that, at least in 
the minimal models with two heavy neutrinos, the strictest bound in the \ie 
channel will put an upper limit on Im(z) which suppresses also the other re 
and Tfj, channels. 

3 Constraints on Yukawa couplings from \l — e 
conversion in nuclei 

The strictest bound on the \ieZ coupling is obtained by the \i — e conversion 
in a nucleus. The current bound on Br^ e = T conversion /T ca p ture was set by the 

(Ti) 

SINDRUM collaboration from the experiments on titanium with Br^e < 
4.3 x 1(T 12 [26] and gold target setting the Br { ^ u) < 7 x 10" 13 [27], both at 
90%CL. 
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Nucleus 


V^[m 5 J 2 ] 




TcaptMreflO S 






0.0396 


0.0468 


2.59 




0.0974 


0.146 


13.07 



Table 2: Data taken from Tables I and VIII of |28j . 



To get the constraint in the fie channel from these experiments, one needs 
to know the expression for the rate in different nuclei. A detailed numerical 
calculation has been carried out by [2H] and we use their formula in Eq.(14) 
to calculate the desired conversion rate. The dominant contribution is due 
to a tree-level exchange of the Z boson, the tree level Higgs amplitude being 
suppressed by the smallness of the charged lepton masses. Other contribu- 
tions, involving also the singlet Yukawa couplings, are suppressed by a loop. 
Therefore at the leading order, the rate depends on the vectorial couplings 
only and using the notation of |28j we have 



can version 



2G% 



where Gf is the SM Fermi coupling and g^'^y are found to be 



~(p) 
9lv 

~(n) 

9lv 



2 (3 - Asl) L 



12- 



= 2 (3 - 44) R 



z 

12) 



QL 



12' 



~(n) 

9rv 



6R 



12- 



(9) 



(10) 
(11) 



Throughout the paper we use L w,z (R w ' z ) to denote the left (right) handed 



couplings of the fermions to the gauge bosons as defined in Eq.(A.15) and 
s w = sind w , where 6 W is the weak mixing angle. The values of depend 
on the given nucleus and are specified in Table I of |28j, while the capture 
rates are given in their Table VIII and we list the relevant quantities in Table 
[2j The resulting bound on the LFV couplings is 



\T Z \ 2 + 

\ L 12\ + 



R 



Z |2 
12 I 



< 10 



-14 



10 



-15 



(12) 



for Ti and Au, respectively. Note that the bound due to the more recently 
measured Au channel is an order of magnitude stronger than the Ti bound 
previously considered in the literature [T6l [18]. After allowing to vary the 
poorly known neutrino mass parameter #13 within the allowed range in table 
1 and the unknown phases S and 0, we obtain in the minimal models a bound 
on lm(z) < 7.5(7.1) for normal (inverted) hierarchy in case of one triplet and 
one singlet and lm(z) < 7.2(6.8) for two triplets, all at the reference mass of 
rriT = 100 GeV for the lightest triplet. 
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When Im(z) is so large, the branching ratios for decays of the triplets to 
light leptons are fixed by the neutrino mixing parameters and can be checked 
at the collider, if the triplet is light enough to be produced [II] . On the other 
hand, the same constraint puts an upper bound on the Yukawa coupling of 
the singlet, which makes it very hard to observe, even if it were light. 

4 Other constraints 

Previous phenomenological analyses of various experimental constraints on 
type III see-saw models [HJ HH [18] considered the three charged lepton flavor 
transitions separately. As explained above, such treatment is not necessarily 
justified, since the relative strengths of the various flavor transitions in type 
III are governed by neutrino mass and mixing parameters - not all entries 
in the Lfj w , Rfj W coupling matrices are independent. Consequently the 
impact of the various constraints should be compared through their bounds 
on the remaining free parameters of the model. In the case of the minimal 
III and I+III models this is the single complex parameter z. As explained in 
the next section, the general observation remains valid even in non-minimal 
type III models, albeit with more free parameters to be constrained. In the 
following we will compare the current and prospective bounds on the Im(z) 
in the minimal III and I+III models. These are to be contrasted with the 
benchmark limits set by \x — e nuclear conversion experiments. In most cases, 
the differences in bounds obtained with normal or inverted hierarchies and 
between III and I+III setups are not significant given their overall size. The 
bounds are not very sensitive to the angle #13 and the Dirac phase 5 due to 
the smallness of #13. There is a mild sensitivity to the Majorana phase </> as 
shown on figure [T] where the comparison of various bounds is summarized for 
the minimal type III case. 

Also to be kept in mind is that the dependence of observables on the 
size of lm(z) is exponential, so that an improvement of a particular bound 
on Im(,2) by 0(1) requires (numerically roughly two) orders of magnitude 
improvement in the actual experimental limit. On the other hand the LFV 
and lepton flavor universal (LFU) effects decouple quadratically with the 
lightest triplet mass as shown explicitly in appendix [A] 

Finally, the situation can also be viewed from the opposite perspective. 
Since within the I+III setup, the r — fi and r — e LFV transitions for example, 
are constrained by the fi — e bound, any positive indication of the other 
transitions in the near future would indicate LFV contributions beyond the 
minimal models. 

Leptonic LFV decays are closely related to the fi — e conversion processes, 
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Figure 1: Comparison of various LFV and LFU bounds on the minimal type 
III model for normal (top) and inverted (bottom) neutrino mass hierarchy. 
The bounds coming from \i — e transitions are plotted in red, r — e in blue and 
T — fi in green. Constraint from the Z width to electrons is shown in magenta, 
to muons in cyan and to taus in yellow. The bounds constrain lm(z) at the 
reference triplets' mass of 100 GeV and depend on the unknown Majorana 
phase 0. Dependence on the other poorly known neutrino parameters is 
negligible as explained in the text. 
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since both receive dominant contributions form tree level Z exchange in type 
III models. The relevant decay widths at leading order are 



q2 

V li-^l^ktA = (\ L ij\ + \ R ij\ ) (\ L kk\ + \ R kk\ ) , ( 13a ) 



+ 2 (\ L ?3 L fi\* + \ R ij R fj I 2 ) 



(13b) 



where we have neglected the final state lepton masses and doubly flavor 
suppressed amplitudes. 

Taken the experimental limits from [21], it turns out that these decays 
constitute the most sensitive bounds on Im(z) coming from r — £ transitions. 

Radiative decays of the charged leptons also put limits on the LFV cou- 
plings. Since the photon coupling to the leptons remains universal at tree 
level, this process has to go through a loop. We have calculated the ampli- 
tudes coming from the W, Z and Higgs loops and we give the result in the 
appendix [B] 

As seen in figure [TJ the limits coming from the loop suppressed \i — *■ 
decay are substantially weaker than fi —>■ 3e, also due to a better experimental 
bound for the latter. 

Semileptonic LFV tau decays r — > tt°£ and r — > T)£, where £ — fi, e were 
identified in [HI [TH] as promising LFV signatures in the tau sector at low 
energies. Present experimental limits on the branching ratios are at the 10 -8 
level [21] and thus the corresponding bounds could be in principle comparable 
to the ones from r — > 3£. The decay widths induced by generic LFV Z 
couplings can be written as 



^2 f2 



(1- ml/ml) 2 (\Lif + \R!f) , (14) 



where h = vr ,^ and fh is the corresponding decay constant. 

The above formula neglects final state lepton masses, but is accurate to a 
percent level even for the r —>■ rjfi channel. In our numerical analysis we use 
the complete kinematic formula, which can be found e.g. in [18] FJ In our 
treatment of the hadronic matrix elements of the rj we follow the formalism 
of [22] and sum over contributions from all light quark flavours (dd, uu and 
ss) as also done in [18J. The derived limits on lm(z) coming from tt° and rj 



1 We could reproduce all the 7r and rj bounds in Table 1 of [TH] except the one for r — > ire 
for which we instead get in their notation |e er | < 6.0 x 10~ 4 . Therefore, this bound is not 
stronger than the one from r — > 3e decay as claimed in the paper. 
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channels differ only slightly and exhibit identical 4> dependence, therefore we 
do not plot them separately in figure [T] 

Z decay widths to lepton pairs of various flavors were measured at LEP 
|24j . Both flavor conserving (diagonal) as well as flavour changing (off- 
diagonal) decay modes could impose relevant constrains on deviations from 
the universal Z couplings to leptons. The relevant leading order decay width 
formula is 

r z _^. = G / m £ ((1 - {m k - m tj ) 2 /m 2 z ) (l - (m A + m tj ) 2 /m 2 z )) [ 
buy 2, 

(\Lf 3 \ 2 + I4I 2 ) (2 - (m| + m|)/m| - (ml - m|) 2 /m|) + (15) 
12 Re (LlRf*)m li m l ./m% 

where the finite lepton mass effects are only important for the tau channels. 

Comparing to experimental measurements listed in [24], presently, the 
flavour diagonal channels yield bounds comparable to those from r decays. 

Charged current lepton flavor universality tests can also probe for signs of 
violations of unitarity of the coupling matrix between light leptons, present 
in I+III models (the 3x3 submatrix of Lfj). The most relevant observables 
here are (semi)leptonic kaon, pion and tau decays, while direct W decay mea- 
surements at LEP yield somewhat weaker constraints. A model-independent 
analysis was performed in [16J. The best bounds on the deviations of uni- 
tarity are at the level of a few per-mille. When translated onto the bounds 
of the I+III model parameters, these are already quite weaker than other 
aforementioned constraints. 

Additional constraints studied in the literature include anomalous mag- 
netic moment of the muon [23], LFV leptonic and semileptonic decays of 
mesons, muonium - anti-muonium oscillations |18] . all of which yield much 
weaker constraints than the ones mentioned above. 



5 Beyond minimal models 

Before concluding, let us comment on non-minimal models with more than 
two heavy fermions. In the minimal cases above, we had 11 real parameters 
governing the Yukawa couplings: two Majorana masses of the heavy fermions, 
5+2 parameters (the PMNS matrix and two masses of light neutrinos) mostly 
fixed from the oscillation data and finally, a single complex angle z which 
specifies a complex orthogonal matrix O. Extending the model with another 
heavy fermion brings in another mass and another phase in the PMNS and 
also a third light neutrino mass (since the overall scale is unknown) and we 
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now have 3 complex angles which specify the 3 by 3 orthogonal matrix O, 
altogether 18 parameters. 

Although there are more free parameters in this case, correlations be- 
tween different channels are generically preserved. This can easily be seen 
by considering the non-universal coupling, 



= ^ y\J m i mU jl m a) O ai O aj U e iU^, 

a=l i,j=l 

where we sum over all the elements of the orthogonal matrix O, regardless 
of the flavor. Therefore one cannot easily enlarge the r£Z couplings by 
enhancing a single element of O without affecting the \xe channel and running 
in contradiction with the fi — e conversion experiment unless one aligns (fine- 
tunes) the available phases. This result holds for an arbitrary number of 
additional triplets and shows that the overall rate of the flavor processes is 
naturally dictated by the most constraining channel. 

On the other hand, there is a potential gain in considering non-minimal 
models with three extra triplets. Namely, one can use the freedom of setting 
the overall scale of neutrinos at will and a positive signal is possible even for 
natural values of the Yukawas. For example, if light neutrinos are degenerate 
with the sum of their masses close to the upper limit from (3 decay and 
cosmology (say < eV [32]), present \i — e conversion bounds already 

probe values of Im(zj) ~ 2 — 4. 



6 Conclusions and outlook 

The fi — e conversion limits will be further improved in the future by several 
orders of magnitude. According to proposals (29] and [50} 151]. one can expect 
a sensitivity of 10~ 16 or even 10~ 18 by the PRISM/PRIME experiment. Such 
a sensitivity would constrain Im(z) to 4.2 (3.8) in case of the minimal I + III 
model and to 3.8 (3.4) for the minimal type III, again for normal (inverted) 
hierarchy. For non-minimal models with degenerate eV scale neutrinos, these 
experiments would already probe Im(zj) < 1. Since the imaginary values of 
Zi are free parameters of the model and setting any of them to zero does not 
enhance the symmetry of the Lagrangian, we consider such values natural. 
We plot both projections in figure [2] against the maximum value of Im^) in 
non-minimal models. 
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0.5 1.0 1.5 2.0 2.5 3.0 3.5 
Max[Im(z,)] 



Figure 2: Present and projected sensitivity of /i — e conversion experiments 
in non-minimal type III see-saw models for a massless lightest neutrino (in 
red empty circles, for both hierarchies) and for degenerate scenario at 1 
eV (in blue filled spades). Minimal model predictions are drawn in green 
dashed lines. In all cases we put all the Majorana phases to zero and vary Zi 
randomly. 

It is worth contrasting this with the projected sensitivity of the MEG 
experiment for the branching ratio of /x — > ej, which is of the order 10~ 13 
|33j . This corresponds to probing values of Im(z) ~ 8.8 in the minimal 
models. 

Another interesting feature of a potential \i — e conversion signal would 
be its ability to distinguish the type III contribution from type I and II. 
This is due to the fact that the dependancy of the vectorial gauge boson 
couplings ^/( p ' n ) on the number of protons in the nucleus is different from 
the contributions of a scalar or a dipole operator (see |28j for details). These 
may be present in type I and II, however they are loop suppressed and we do 
not consider them here. By measuring the conversion rate of two different 
nuclei, one can check for the vectorial nature which should dominate for the 
type III contribution. 

The tau LFV decays are expected to be improved at a future flavor factory 
by one to two orders of magnitude [35J . The Z leptonic width measurements 
could possibly be improved at the LHC or a future linear collider although 
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we are not aware of existing dedicated sensitivity studies. On the other hand 
LFU tests in charged currents will be difficult to improve due to limiting 
theoretical uncertainties, although BES111 could improve on the present ex- 
perimental precision [34] . 
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A Derivation of the fermion couplings 

In this appendix we derive general expressions for the couplings of the charged 
and neutral fermions to the SM gauge and Higgs fields in the presence of a 
fermionic singlet and a weak triplet with hypercharge 0. We start with the 
Lagrangian in a two component notation and derive the rules in the four- 
component notation in the physical mass basis. In order not to clutter the 
notation, we initially consider an addition of a single fermionic triplet and a 
single singlet and we generalize the result for an arbitrary number of triplets 
{tit) and singlets {ris) in the end. 

The starting point is the Lagrangian written with two component Weyl 
spinors in a basis where the Yukawa matrix of the charged fermions is real 
and diagonal and the Majorana masses ttit and mj are also real 



Here, i, j are the family indices running from 1 to 3 and the standard covari- 
ant derivatives are defined as 



C e = i 



y'/rfLitj + y l s H T ia 2 LiS + y l T H T ia 2 a a T a Li + h.c. 
1/2 (m T T a T a + m s SS) + h.c. . 



(A.l) 



D L » = d,- tg/2AX - ig'/2B, 
D T » = d, + ige abc T h A^ 



(A.2) 
(A.3) 
(A.4) 
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where a, b, c are the usual SU(2) indices and cr's are the Pauli matrices. 

We use the linear combinations of the fields, labeled by their U(l) charge 
\plT ± = T 1 =p iT 2 and T° = T 3 and after spontaneous symmetry breaking, 
the Higgs field becomes 



H 



(v + h + ix)/V2 



(A.5) 



and the Lagrangian in Eq.(|A.l|) gives the following mass terms 
r 

f — f m,as 



-@ T+)M,(^_)-( 



T S)M U [ T | /2 + h.c. 

S 



(A.6) 



where 



' y ) and M v 



v y T 



rrix 



3 x3 v Vt v y l s > 
' >•• 



K v y 3 s 



(A.7) 



m 5 



can be brought to a diagonal form by a biunitary and congruent transforma- 
tion for the charged and neutral fields 



M, = U +] M e U- 



U 0T M U U° . 



(A.8) 



In the limit when vtjt "C m?, one can expand these matrices in terms of 
small parameters £j = vy^/rriT, £si — v y l sl m Si an d s\ = v y^rrii/m 2 ^, 



u° 



A-IKI 2 o 







1-Me'\ 2 












2 p l * 


\ -e' -e' 


-4 


l-T- 


/I _ Ik |2 _I £ *£ 


2 c e c r 


e* 


_l e 1 _ I| e |2 


-Ve 
2 v T 


e* 


— -F F* —-F F* 
2 c e c T 9 M T 


l-\\e T 


2 * 








/ _ + e* Si e S k))U k j 





\ 



(A.9) 



(A.10) 



—£kUkj/ V2 
-£skUkj/V2 



After performing these rotations, we combine the mass eigenstates of the 
charged fermions and the triplets into a four component Dirac spinor while 
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the neutral fermions form a Majorana spinor using the usual prescription 

^- (£).*- ft). *■-(£)■ (-) 

The mixing matrices alter the gauge couplings of the SM fermions and since 
they mix the chiral fermions with vector-like triplets, it is convenient to 
introduce a general notation for the charged and neutral four component 
spinors (see also the appendix of [14J) 

/- = (e,/x,r,T-), $ = (is 1 ,v 2 ,v 3 ,T* i ,S). (A.13) 

Using such a convention, we can write down the W and Z couplings in a 
unified way with 

y = (£ ;')/=ftV),,' = (; _° 1 ),ft,H = ^, (A.i4) 

and we have a Lagrangian 

£ int = -e/^/i + (gj'^ + {L w P L + R w P R ) ij f j + h.c.) + (A.15) 

(0 + 7; (L*P L + R*Pr) 31 fi + h.c.) + (A.16) 

f /^(L Z P L + J^P*)^ + X/, (^ x Pl + P*Pi4- /i (A.17) 

+ hf l (L h P L + R h P R ) tj f J , (A. 18) 
with the following gauge 

^ = U%l%/y/2 + U%Up p R* = U%U+, (A.19) 

(A.20) 

and would-be-Goldstone and physical Higgs couplings 

L* = y?- 3a (yfil&lfe - U%U^) + yT^K^ (A. 21) 

K = -y?u° a *u+, (A.22) 

L% = ^V?U+*U^ - iyZ-^UffUr,, (A.23) 

l _ v <*TT-*TT+ 1 ^,P~ 3a *l 
1 

1 



Rl = -^y?U-*U+ + iyt 3a *U-*U+. (A.24) 
4 = -^y?u+*u^ - vfr^UffU^, (A.25) 

4 = -\yWatK ~ yT~ 3a *u-*u+. (A.26) 
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In the notation above, repeated indices are always summed over. The indices 
a, a' run over the light families from 1 to 3, /3 runs over the number of 
triplets from 4 to 3 + tit, while 7 is the singlet index going from 4 + nj> to 
3 + riT + ns- When additional copies of particles are considered, the mass 
matrices in Eq.(A.7) have to be extended. 

A couple of features of the model are noteworthy. While the photon 
vertex remains universal at the tree-level, the Z vertex now receives off- 
diagonal entries. Also, the right-handed couplings are now present, however 
they are always suppressed by the mass of the light charged fermions m^/rriT 
which can be seen from the expansion of U + in Eq.(A.9). Notice that the 
SM limits are easily obtained, by either sending yx,s — > and/or m^s — ► 00. 
In this case, the mixing matrices become diagonal and the SM expressions 
are recovered. 



B ji —> /17 calculation 

Here, we discuss the calculation of amplitudes for the f 2 — > /17 decay. We 
have done the calculation in gauge with arbitrary left and right-handed 
gauge couplings of the fermions with arbitrary masses m 12 of / 1)2 . The 
amplitude is proportional to the d = 5 operator 

i< V e'V, (B.l) 

where e is the polarization vector of the photon with momentum q = p 2 — p\ 
and Pi are the four-momenta of /j. The final result has to be finite and £ WjZ 
independent. 

Before giving the transition amplitude, we would like to comment on the 
divergency cancellations in models with non-unitary mixing matrices which 
is the case for the type III seesaw. If the W coupling matrix is unitary, the 
divergent part proportional to ^ vanishes or it is cancelled by the diagrams 
with photons radiating from the external fermions. This does not happen 
in models where vector-like fermions mix with the chiral. The problem is 
resolved by noting that 

{L W *L W ) .. = -L* i^j, (B.2) 
which is a consequence of the SU(2) structure of the electroweak Lagrangian. 



The relation in Eq.(B.2) holds also for R couplings, and both can be checked 



from Eqs.(A.19) and (A.20). Given that the non-unitarity of the mixing 
matrix is directly related to the non-universality of the Z coupling, we expect 
the cancellation to come from a diagram with a single off-diagonal Z coupling. 
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Indeed, when we calculate the Z — 7 mixing diagrams, the divergent part 
vanishes. 

In order to get a finite and gauge invariant result, we have to sum the 
diagrams with unphysical would-be-Goldstone fields. Their couplings can be 
related to the gauge boson couplings (see also [36]) 

L tj = g/mw { L Yj m i - R Vj m i) , ( B -3) 
Lfj = g/ic w m z (Lfjirii - RfjUij) , (B.4) 

and symmetrically for R — > L. Using these relations, all the £ dependent 
terms cancel to all orders in 7711,2 and we have a finite, gauge invariant result 
coming from the W, Z and Higgs loops, together with corresponding (p and 
X loops. We expand the scalar integrals in small m^-, set m\ = and get the 
amplitudes 



q 3+n T +n s -| 

v n=l 



Ml 



M 



M' 



Qm n L%*R% (xl - I2x 2 n + 6a£ log x n + 15x n - 4) (1 - x n ) + 
m 2 L%*LW (4a£ + 18x 3 n ^gx n - 49^ + 78a£ _ 43x n + 10) 

3+n T 

v c=l v 

6m c Lf c R^ 2 (xl — 6x c \ogx c + 3x c — 4) (1 — x c ) + 
m 2 Lf c L z c2 (bx 4 c - Ux\ - l%x 2 c \ogx c + 39x^ - 38x c + 8) 

f 1 ia ftv e^P R h J2 77 

c=l ^ 



967r 2 m^ 



Qm c R h lc R h c2 (y 2 -Ay c + 2 log y c + 3) (1 - y c ) - 
m 2 R h lc L h c2 (yl - 6y 2 c + 6y c log y c + 3y c + 2)) 



(B.5) 



(B.6) 



(B.7) 



where n sums over the neutral particles in the loop (three light neutrinos 
and riT + n$ heavy mediators), c sums over the charged particles (three light 
e, /1, r and n T heavy triplets) and x n = m^/rriy^, x c = m 2 /m 2 z , y c = m 2 /m\. 
The amplitude proportional to Pl is obtained by substituting (L,R) — > 
(R, L). Finally, the total decay rate for the process is given by 



16/T 



\M L \ 2 + \M R \ 2 ), 



(B.8) 
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with M l ,r = Ei=w,z,h Mi tR . 

Notice that the result above is valid also for theories with right-handed 
gauge couplings, e.g. in left-right symmetric theories. From the results above, 
one can easily reproduce the calculations for type I case. The Z and H 
amplitudes are zero, so are the R couplings, therefore the only contributing 
piece is the third line of Eq.flR5]), proportional to L-J = U$/y/2. With this 



substitution we reproduce the well-known results in [37]. We cannot fully 
reproduce the results of [3E] for the case of pure type III, our result for the 
decay rate is bigger by roughly a factor of two. 
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